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Comment on astrophysical consequences of a neutrinophilic 2HDM 
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Several authors have pointed out that the scalar-mediated interaction of neutrinos in a neutro- 
philic two-Higgs-doublet model (V2HDM) can be as strong as electromagnetic interaction 0-0]. 
We show that the coupling constants of neutrino-scalar interaction are actually restricted to be 
Hi ^ 1-5 x 10 -3 by supernova neutrino observation, and further constrained to be yt < 2.3 x 10 -4 
by precision measurements of acoustic peaks of the cosmic microwave background. Based on the 
energy- loss argument for supernova cores, we derive a slightly more restrictive bound yi < 3.5 x 10 -5 . 
Therefore, the i^2HDM has lost its spirit of explaining tiny Dirac neutrino masses while keeping neu- 
trino Yukawa couplings of order one. 

PACS numbers: 12.60.Fr, 14.60.Lm 



The z^2HDM extends the standard model with three 
right-handed neutrinos and introduces exclusively for 
neutrinos an extra Higgs doublet tp, which acquires a 
small vacuum expectation value (vev) v v =0.1 eV. In 
this model, neutrinos are Dirac particles, and their mass 
matrix is given by M v = Y v v v with Y v being neutrino 
Yukawa coupling matrix. A salient feature of the Z/2HDM 
is that Yv can be of order one even for sub-eV neutrino 
masses [JH2I . In the flavor basis where the charged- lepton 
Yukawa coupling matrix is diagonal, one can further take 
Y u to be Hermitian by rotating right-handed neutrino 
fields in the flavor space. Hence the Yukawa interaction 
of neutrinos can be written as 
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where V^Y U V = Diag{j/!, y 2 , 2/3} with V being the neu- 
trino mixing matrix, which relates neutrino mass eigen- 
states Vi to flavor eigenstates v a , and rrij = yiV v (for 
i = 1,2,3) are neutrino masses. Here 77 is a scalar boson 
arising from the neutral component of tp, and its mass is 
naturally around the vev of tp, i.e., m v s» v v = 0.1 eV. 
Some cosmological and astrophysical consequences of the 
neutrino interaction in Eq. (1) with yi ~ 0(1) have been 
discussed in Refs. [ij-yl], however, the restrictive bounds 
on neutrino Yukawa couplings are unfortunately missed. 

In fact, stringent bounds on the neutrino-Majoron in- 
teraction have been obtained in the literature by assum- 
ing a pseudoscalar coupling iyiyl^ViX with \ being a 
pseudoscalar boson 0, Q. One can show that those 
bounds apply as well to the scalar case in the relativis- 
tic limit, where small neutrino masses can be neglected. 
However, it should be noticed that the lepton-number- 
violating processes are forbidden in the v2HDM. 

The observation of neutrinos from Supernova 1987A 
requires that the mean free path of electron antineutrinos 
in the presence of cosmic background particles should be 



larger than the supernova distance, i.e., D < 1 with 
D = 51.4 kpc, in order to avoid significant reduction of 
neutrino flux Q. The relevant processes are v e + 77 — > 
v e + r/, v e + v e — > rj + 77 and v e + v a v e + v a for 
a = e, t. After removing the lepton-number-violating 
contributions from the Majoron model [J], one can obtain 
a restrictive bound on neutrino Yukawa couplings 



V, 



< 1.5 x 10" 3 



(2) 
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If 77 bosons decay rapidly into neutrinos and are absent 
in the cosmic background, the bound will be weaker but 
on the same order of magnitude. Given v v — 0.1 eV, 
neutrino Yukawa couplings are determined by neutrino 
masses yi = rrii/v^, thus heavier neutrino mass eigen- 
states interact more strongly with the scalar boson. In 
the case of normal mass hierarchy, the bound in Eq. (2) 
may be slightly relaxed to j/j < 10 -2 , because electron 
antineutrino possesses a small fraction of the heaviest 
mass eigenstate@. For the inverted mass hierarchy or 
nearly degenerate neutrino mass spectrum, however, the 
bound in Eq. (2) is still applicable. 

As indicated by precision measurements of the acous- 
tic peaks of the cosmic microwave background, neutri- 
nos should be freely streaming around the time of pho- 
ton decoupling T T)( jec = 0.256 eV in order to avoid 
the acoustic oscillations of the neutrino-scalar fluid [f|. 
At this moment, the neutrino temperature is T^dcc = 
(4/11)1/3^ doc _ 0.183 eV. For the relevant two-body 
scattering processes v\ + r/ — > z/j + 77, z/j + — > 77 + 77 
and Vi + Vi — > Vj + 9j via 77-exchange, we can sim- 
ply estimate the scattering rate as T u s» yfT U: dec U P to 
some numerical factors. The free-streaming argument 
requires T„ to be smaller than the cosmic expansion 
rate H lAcc = 100 km s" 1 Mpc' 1 (fluh 2 ) 1 / 2 (z dcc + l) 3 / 2 
with Vlyih 2 = 0.134 being the cosmic matter density and 
Zdcc = 1088 the redshift at photon decoupling. Hence 
one can derive a more restrictive bound < 1.1 x 10~ 7 
5]. Taking account of the existing bound on Yukawa 
couplings in Eq. (2), one should increase v v by three or- 
ders of magnitude to guarantee sub-eV neutrino masses. 
Therefore, the mass of 77 is expected to be « v v = 
100 eV, which is much larger than neutrino temperature 



2 



at the time of photon decoupling. As a consequence, 
77 bosons have already decayed into neutrinos, and the 
relevant process Vi + vi — > Vj + Dj is mediated by a vir- 
tual 77 boson Q. The scattering rate is modified to be 

r ^ ~ yt T u,deo/ m t ~ Vi T u,dec/ V t = vt^dcc/ ' m t '> thuS 

the true bound from the free-streaming argument is 

Vl < 2.3 x 1(T 4 , (3) 

for rrii ~ 0.1 eV. Since the scattering rate is propor- 
tional to yf , the neglected numerical factors are indeed 
unimportant for the bound in Eq. (3). 

An important point is that 77 is massive enough to de- 
cay into neutrino- antineutrino pairs 77 — > Vi + Pj. The 
lifetime in its rest frame is t v = (3yf m^/l&ir)^ 1 w 
1.1 x 10~ 9 s, where yi — 10~ 4 and m v ~ v v = 1 keV have 
been taken. Although rj bosons can be copiously pro- 
duced in the supernova core, one may expect that they 
will decay soon after production and thus cannot cause 
excessive energy losses. However, since the temperature 
in the cooling phase is sufficiently high T = 30 MeV, the 
lifetime of thermal 77 bosons should be lengthened by a 
Lorentz factor E/m v w 10 5 . Consequently, the relativis- 
tic 77 bosons before decaying may have traveled a distance 
l v R3 3.3 x 10 6 cm, which is larger than the core radius 
R = 10 km. But 77 bosons cannot freely propagate in the 
neutrino background, their mean free path can be esti- 
mated as A „ = (yfT)- 1 w 6.6xl0 3 cm that is comparable 
to the mean free path of neutrinos in the case of standard 
neutral-current interaction. Thus 77 bosons behave like a 
new species of neutrinos and accelerate the energy trans- 
fer, which leads to the reduction of the cooling time or 
the duration of supernova neutrino burst. There are two 
possibilities to avoid the contradiction with the neutrino 
observation of Supernova 1987A [3]: (i) to increase the 
coupling (yi > 10 -4 ) such that A,, becomes much smaller 
and the energy transfer by 77 bosons is negligible; (ii) to 
decrease the coupling (yi < 10~ 4 ) such that 77 bosons 
have never been trapped and thcrmalized in the core. 
The first possibility is already excluded by the bound in 
Eq. (3), while the second one is subject to the constraint 
from standard energy-loss arguments. 

The production of 77 bosons will be efficient via the 
bremsstrahlung process Vi + N — s- f,- + N + 77 because 
of the high nucleon density. The emission rate should 
be proportional to the thermal average of vN scattering 
rate an v n B where a is the vN collision cross section, n v 
and 71b are respectively the neutrino and baryon number 



densities. The energy of emitted 77 bosons is of order T. 
Put all together, we can get the volume emission rate 

Qn » 54yfGlT 6 n B « 2.4t/, 2 x 10 42 erg cm- 3 s" 1 , (4) 

where a = G\E 2 , T = 30 MeV and n B = p/m N with 
p = 3.0 x 10 14 g cm -3 have been assumed. Note that 
all the neutrinos have been taken to be relativistic and 
non-degenerate, which is an excellent approximation for 
and v T . For degenerate v e , there will be a blocking 
factor that suppresses the scattering rate, which has been 
neglected in Eq. (4) for an order-of-magnitude estimate. 
The energy loss should be small so as not to shorten the 
neutrino burst, so we require the volume emission rate to 
be Q v ^ 3.0 x 10 33 erg cm~ 3 s~ 4 and then obtain 

i H < 3.5 x 10- 5 . (5) 

As a matter of fact, there are additional contributions to 
the production of rj bosons, such as z/j + e~ — > vi + e~ + 77 
and Vi + Vi — > r\ + 77. Hence the bound may be slightly 
stronger if all the contributions are included. One can 
then verify that the mean free path A^ ~ 10 s cm is much 
larger than the core radius, which is consistent with the 
prerequisite that 77 bosons can escape from the core and 
carry away energies. 

Based on the above discussions, one may take yi = 
10 -5 and figure out the cross section of Vi + i>i — > Vj + Pj 
via 77-exchange, and that via Z°-exchange. It is straight- 
forward to get a v ~ yf/E 2 for the former case, while 
<?z ~ GpE 2 for the latter, where E is the neutrino en- 
ergy. In the neutrinosphere with a typical temperature 
T = 10 MeV, we further obtain o^/oz « yf/G^E 4 w 
10~ 4 for Vi = IO- 5 and E = 3T = 30 MeV. Therefore, 
the scalar-mediated neutrino interaction is too weak to 
equilibrate supernova neutrinos of different species in the 
neutrinosphere and thus cannot wash out the collective 
effects in supernova neutrino oscillations [l[ . 

In conclusion, if the astrophysical bounds on the scalar- 
mediated neutrino interaction are taken into account, 
the motivation for the i^2HDM becomes very weak in 
the sense that tiny Dirac neutrino masses cannot be ex- 
plained by a small vev but large Yukawa couplings. 
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